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Abstract. Polar magnetospheric gaps consume a fraction of the electric potential that develops accross open field lines. 
This efi"ect modifies significantly the structure of the axisymmetric pulsar magnetosphere. We present numerical steady-state 
solutions for various values of the gap potential. We show that a charge starved magnetosphere contains significantly less electric 
current than one with freely available electric charges. As a result, electromagnetic neutron star braking becomes inefficient. We 
argue that the magnetosphere may spontaneously rearrange itself to a lower energy configuration through a dramatic release of 
electromagnetic field energy and magnetic flux. Our results might be relevant in understanding the recent December 27, 2004 
burst observed in SGR 1806-20. 
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1. Introduction 

The magnetosphere of a rotating neutron star with polar mag- 
netic field B„ mass M, ~ Mq, radius r, ~ 10km, magnetic 
dipole moment fi - rlBt/2, and angular velocity Q. is expected 
to radiate electromagnetic energy at a rate 

iem ~ a . (1) 



a is a factor of order unity (IBeskin 1997l l'. The source of the 
radiation is the neutron star rotational kinetic energy which is 
lost at a rate 



4 2 . 

i'kinetic ~ -M.r.Qfi . 



(2) 



It is usual to equate eqs. ^ & |2] and thus obtain an estimate 
of the stellar magnetic field B,. In general, however, the two 
do not have to be equal. As we shall see below, in the case of 
axisymmetry, electromagnetic torques need to be significantly 
revised. 

In the context of ideal axisymmetric MHD, electric charges 
are available in copious amounts and move freely along mag- 
netic field lines, shorting out any component of the electric field 

* Present address. -Reseaich Center for Astronomy 
& Applied Mathematics, Academy of Athens, e-mail: 
icontopOacademyofathens . gr 

' a = ^ sin^ 6 for a misalligned dipole rotating in vacuum. In 
that scenario, an alligned magnetic rotator (d ~ 0) does not radi- 
ate. However, when the neutron star is not surrounded by vacuum, 
one needs to consider the structure of its rotating charged relativistic 
Goldreich-Julian-type magnetosphere (Goldreich & Julian 1969). In 
that case, the electric currents that flow through the magnetosphere 
lead to electromagnetic energy losses comparable to the ones for a 
misalligned magnetic rotator. See the Appendix for a general calcula- 
tion. 



that might arise along the magnetic field. As a result, magnetic 
flux surfaces become equipotentials, and an electric field E de- 
velops accross magnetic field lines (E ■ B = 0) with magnitude 



-Bn 



(3) 



where, Q.f, a constant along magnetic flux surfaces (see below), 
can be thought of as the angular velocity of rotation of mag- 
netic field lines (r is the cylindrical radius; Bp is the poloidal 
component of the magnetic field). The source of the electric 
potential accross magnetic field lines is obviously the rotating 
magnetized conducting surface of the neutron star which acts 
as a unipolar inductor. For the sake of simplicity, most stud- 
ies of the axisymmetric pulsar magnetosphere have assumed 
that the full potential drop induced accross field lines along the 
surface of the star continues to manifest itself all along those 
field lines, i.e. = Q. It has been pointed out, however, (e.g. 
Ruderman & Sutherland, 1975) that 'open magnetic field lines 
play a role analogous to that of conducting wires in ordinary 
circuits. If the wire is broken near the pulsar surface, a po- 
tential drop develops accross the gap'. The presence of such 
gaps obviously reduces the electric potential available accross 
open field lines, and thus the electromagnetic energy power ra- 
diated at large distances. Models of particle acceleration and 
pair creation above the polar cap of rotation-powered pulsars 
yield potential drops near the surface of the star of the order 
of 10'^ Volts (e.g. Hibschman & Aron s 200U Arons personal 
communication), and therefore, in general, Q/r < Q. 

Beskin & Malyshkin 1998 took the above well known ef- 
fect into account in their calculation of the modified magneto- 
spheric structure inside the light cylinder In the present paper 
we obtain the first global solution of this problem. In § 2 we 
formulate the problem and the numerical method that we im- 
plement for its solution. In § 3 we obtain the structure of the 
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magnetosphere for various values of Q/r in the range [0, Q] and 
argue that the magnetosphere may switch between solutions, 
releasing energy in the process. In § 4 we discuss the relevance 
of our results in understanding the recent December 27, 2004 
SGR- 1806-20 burst. Our conclusions are summarized in § 5. 

2. The differentially rotating magnetosphere 

We will work in cylindrical spatial coordinates r, (p, z, and will 
consider only the axisymmetric case where the magnetic dipole 
axis is aligned with the axis of rotation. This simplification al- 
lows us to introduce the magnetic flux function i// (tfr/r is the 
^-component of the magnetic vector potential), the poloidal 
electric current function A = A(iff) (the poloidal electric cur- 
rent contained within the magnetic flux surface tfr is equal to 
Ac/2; - A/r), and the magnetic field line 'rotational veloc- 
ity' Q/r = Q.f{i//). The various magnetospheric physical quanti- 
ties are obtained as follows: 



B = -(-(A,,A,l/r,) 
r 



E = -Vtf/^ -(iffr,0,tf/,) , 

c c 

J = -^V X B = --^ (a>,, i/r„. - ^ + -A'lfr,] , 



An 



4nc 



(4) 



(5) 



(6) 



(7) 



Here, and in what follows, t//_x = difr/dx. Also, (...)' = 
d(. . ■)/dt//. When we neglect inertia, force balance requires that 

-J X B + p,E = . (8) 

c 

Following Gruzinov 2005, we take 

c = ^ = Q = 1 , (9) 
and thus eq.|Slbecomes 

(^iff,, + — + lAzz)- — = -AA' + r^Q.fQ.'p{Vifff (10) 

This is a more general form of the pulsar equation than the one 
considered in Contopoulos, Kazanas & Fendt (hereafter CKF) 
where Q.f = Q. = 1 everywhere. 

is related to the magnetospheric potential drop Vf be- 
tween the axis and any magnetic flux surface i// (eq.|5}, namely 

VfW^- TofdiA (11) 
c Jo 

(in units B,rjQ^/2c^). This is in general different from the stel- 
lar potential drop between the pole and the footpoint on the 
surface of the star of the magnetic flux surface tp, namely 



y,((A) = (A/c . 

The difference 

1 r''' 

V(iff) = V,-VF ^ - (l-nf)di/r. 

c Jo 



(12) 



(13) 



is just the particle acceleration gap potential which devolops 
along the magnetic field near the footpoint of the magnetic flux 
surface (e.g. Beski n 19971 1. In the region of closed field lines 
(hereafter the 'dead zone'), there is no particle flow, and there- 
fore there is no need for the formation of particle acceleration 
gaps. We can thus express 



O.F0 ^ 1 along open field lines i// < i/'open 
1 in the 'dead zone' 



(14) 



Q.F0W is determined by the particle acceleration gap micro- 
physics and by the supply of charge carriers at the base of the 
magnetosphere (see discussion in § 4). In our present discus- 
sion, Q.F0W is essentially a free function. For the sake of sim- 
plicity, we take Q.fo - const., as in |Beskin & Malyshkin 199'8| 
The magnetospheric gap potential is, therefore, given by 



y((A) = 



(1 - D.Fo)4' along open field lines ip < ip open 
in the 'dead zone' 



(15) 



V is minimal at the center of the polar cap and increases oc p- 
as we move away from the axis. 

Observational manifestation of the differential magneto- 
spheric rotation is thought to be found in the sub-pulse slow 
drifts accross the pulse profile in almost aligned pulsars (e.g. 
[Rankin & Wright 2003^ . Interpretation of such drifts remains 
still rather sketchy. We speculate that the sub-pulses are asso- 
ciated with the above mentioned magnetospheric gaps present 
around the separatrix between open and closed field lines where 
the need for electric charge carriers is the greatest (as we 
discuss below, this is where closes the electric circuit of the 
poloidal electric current that flows through the polar cap). In 
most cases with observed sub-pulse drifts (~ 1 sec period pul- 
sars) Qfo is expected to be much smaller than Q.. These gaps 
are probably carried around the axis of rotation by the 'friction' 
between the differentially rotating open and closed line regions, 
and thus their observed angular velocity is found to be close to 
Q.. 

As we mentioned before, solutions of the pulsar equation 
exist only for the simplest case Qf = 1 and - 0. Even 
in that case, though, a strong mathematival singularity, the so 
called 'light cylinder' 



He = 1 



(16) 



makes the problem non-trivial. Note that this is just the force- 
free Alfven surface, and only very recently has a numerical 
method been presented which allowed us to obtain a 'smooth' 
solution that fills all space (CKF). The main features of that 
numerical solution (further refined in Gruzinov 2005) are: 

1. The region of open field lines, the so called 'polar cap', is 
slightly larger than the region of static dipolar field lines 
which cross the equator beyond the distance r - I, namely 

'Aopen = 1.23 (17) 

(present calculation^). We remind the reader that i/'open - 1 
for a static dipole, and therefore we see that rotation 'pulls' 
dipolar magnetic field lines out. 



- In CKF, with a much lower numerical resolution, we obtained a 
value of 1.36. Gruzinov 2005 obtained a value of 1.27 with a numer- 
ical resolution comparable to our present one. 
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2. The distribution of poloidal electric current along the open 
field lines is very close to the one in the Mi chel 1974 rela- 
tivistic split monopole solution, namely 



iCKF - 
j ^Michel 

lo 



= -lA (2 - 7^ I along open field lines 

\ '/'open / 



in the 'dead zone' 



(18) 



The electric current distribution deviates slightly from the 
above near if/ ~ i/'open where field lines follow the singular 
shape of the separatrix between the open and closed line 
regions (see fig.|6lbelow). 

3. The return current of the above current distribution, 
A(i/'open) ~ lAopen, flows along the separatrix. This implies 
the presence of magnetic and electric field discontinuities 
accross the separatrix. 

4. In general, the equatorial extent of the 'dead zone' may 
be taken as a free parameter (see section 5, Appendix). It 
is very natural, however, to assume that the 'dead zone' ex- 
tends all the way to the light cylinder^. 

5. Open field lines become monopole-like around and beyond 
the light cylinder. 

6. |B| > |E| everywhere^. 

We are now ready to address the physically more interest- 
ing case Q.'p + 0, in the simplest possible case where Q.p = 
fjfo = const. < 1 in the open field line region, and Q.^ = 1 in 
the dead zone (eq. ll4> . When \p < i/'open, we can rewrite eg. 1101 
in the new spatial coordinates x = Q.Fof and y = fifoZ, 



(1 



AA' 



(19) 



Eq. ^]is the same as our original equation in CKF. We thus 
expect that solutions of eq. [TO] will be very similar to the ones 
obtained in CKF. We would like to emphasize the following 
interesting features: 

1 . As in CKF, it is natural to assume that the corrotating 'dead 
zone' extends all the way to the light cylinder distance, i.e. 
'"c = He = 1 ■ The real mathematical singularity, however, 
is not at the light cylinder, but at a certain distance outside, 
the 'open field light cylinder' 



'"oflc 



n-J > 1 



(20) 



This is where we will apply the numerical iteration routine 
developed in CKF. 



^ Gruzinov 2005 shows that this solution requires infinite magnetic 
fields at the point r = l,z = (in the limit of infinitesimal grid size). 
[Uzdensky 2003 1 and [Lyubarskii 1990| argue against infinite fields and 
thus conclude that the dead zone should end at some finite distance 
inside the light cylinder. 

* This observation counteracts criticism that the assumptions 
of ideal MHD may break down beyond the light cylinder 
jOgura & Kojima 2003||Spitkov sky 2004i. We believe that the source 
of the opposite result presented in Ogura & Kojima 2003 (their fig. 5) 
is due to their numerical boundary condition, eq. 3.3 & fig. 1, namely 
that field lines become horizontal at large radial distances. 



2. We also expect i^open ~ 1 as in previous solutions. 

3. As in CKF, we expect to encounter similar magnetic and 
electric field discontinuities accross the separatrix between 
open and closed field lines. 

4. The r.h.s. of eq.^]is obtained through a numerical iteration 
along the open field light cylinder that guarantees smooth 
crossing of the singularity. Based on our experience, we 
expect this function to be very close to -AA^j^. Therefore, 
to a good approximation. 



A - QfoAcKF : 



(21) 



i.e. A oc Q.pg. Obviously, as Q.fo — » 0, A — > 0. As we 
will see, this result has very interesting implications for the 
electromagnetic torques on the surface of the neutron star. 

Eq. ^1 is elliptic with mixed boundary conditions inside and 
outside the open field light cylinder r = rotic : 

1 . i/r = along r -Q, and i// - i/'open along the equator beyond 
r - 1 (Dirichlet boundary conditions)^. 

2. - Q (i.e. Br - 0) along the equator in the closed line 
region r < 1 (Newman boundary condition). 

3. \jjr - AA' / (2Q.fo) along the open field light cylinder r - 
'"oflc (Newman boundary condition). 

4. Finally, as in CKF, boundary conditions at infinity are 
irrelevant as long as we rescale our spatial coordinates 
to new ones that map the full (r,z) - ([0, oo], [0, oo]) 
space to our finite grid size (rnew.Znew) - ([0, 2], [0, 1]). 
Note that this is not the case for other numerical schemes 
where the integration is constrained within finite spa- 
tial extent ( jOgura & Kojima 2b03| Goodwin et al. 2004; 
Gruzinov 2005). 

The above show that the problem is well defined inside and 
outside the open field light cylinder, and therefore one can ob- 
tain solutions for a general current distribution A = A(iff). The 
two problems are, however, independent, and in general the so- 
lution will be discontinuous at the open field light cylinder, un- 
less one chooses the one poloidal electric current distribution 
A = A(i/r) that will guarantee i//{r = r^ft^,z) = i//{r = r^ft^,z). 
Continuity will also result in the smoothness of the solution 
(see above boundary condition # 3). A{ifr) is obtained as de- 
scribed in CKF by itteratively correcting to a new function 

AA' L = i [^{r = r„fl,,z) + = C^,z)^ 

\ / new 

= i (AA'(^(r = r^-g^, z))oid + AA'(,p(r = z))oid) (22) 

for all grid points along the open field light cylinder In the 
present work the relaxation inside each grid proceeds together 
with the itteration along the open field light cylinder. This im- 
provement over the CKF method allowed for a much greater 



^ As is shown in the Appendix we are in general allowed to arbi- 
trarily choose the equatorial extent of the closed line region. In that 
case, i/'opcn is obtained as a solution of eg. llOl inside the open field light 
cylinder, i.e. it is not an extra free parameter (see Goodwin et al. 2004 
for a different point of view). 




r 



Fig. 1. Q.F - 1 everywhere (CKF; Gruzinov 2005). Thin lines 
correspond to i// intervals of 0.1. i/' = along the axis. The 
dotted line shows the separatrix ip = i^open - 1.23. The mathe- 
matical singularity is at ric = roOc = 1 • 




r 



Fig. 2. £lfo - 0.8 in the open line region, i/^open = 1.23. roflc = 
1.25 
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r 



Fig. 3. Ofo = 0.6 in the open line region. Thin lines corre- 
spond to i/r intervals of 0.1. = along the axis. The dotted 
line shows the separatrix i// - i/'open = 1 -20. The mathematical 
singularity is at roHc - 1.67. 




r 



Fig. 4. flfo = in the open line region. (Aopen = 1-20. roflc = 
i.e. the mathematical singularity is absent in this limiting case. 



numerical resolution and a much faster speed of numerical con- 
vergence! Our numerical scheme consists of an elliptic solver 
with Chebyshev acceleration JPress et al. 19 88 1 over two lOOx 
100 numerical grids joined along the open field light cylinder. 
The discontinuities of A{i//) and Qf(i/') accross the separatrix 
between the open and closed line regions are smoothed out nu- 
merically over a distance dtfr - 0.05 inside the dead zone. 

3. Steady-State magnetospheric solutions 

The various types of solutions of eq. ^| are shown in Figs ^ 
Eland Fig. IB. II Thin lines correspond to i// intervals of 0.1. 
\jj — along the axis. The dotted line represents the separatrix 
'/'open- Fig. 1 shows the CKF case Q.f - 1. Here, i/'open - 1-23 
within the accuracy of our numerical simulation. Figs.|2]&|3] 
show intermediate cases with Q.fo - 0.8 & 0.6 in the open line 
region respectively, i/'open - 1.23 <& 1.20 respectively. Fig. 0] 



shows the limiting case with Q.fo - 0. In that case there is no 
light cylinder singularity that would yield the poloidal electric 
current distribution A(^ir). However, we showed previously that 
the poloidal electric current disappears, since it is obtained as a 
limit of solutions with Q^o — > in the open line region. Here, 

i/'open = 1-20. 

The various magnetospheres show a similar poloidal mag- 
netic field distribution. This result is understood since B, is 
approximately oc r^^ in the equatorial dipole-like closed line 
region, and therefore an approximate estimate for i/'open is 

i/'open ~ — = 1 . (23) 
rc 

However, they differ significantly in the amount of electric 
charge and electric current they contain in the open field line 
region, since pe oc Qfg and A oc Q^^. As a result, they diff'er 
in the amount of electromagnetic field energy they contain in 
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0.4 0.6 
Omega Fo 

Fig. 5. Summary of our numerical solutions applied in the case 
of SGR 1806-20. We show here (Aopen (continuous line), the 
accelerating potential Vacc/ClO'^ Volts) (dashed line), and the 
spindown rate IPI/IO " Hz sec ' (short dashed line). On the 
plot are shown our estimates for the magnetospheric configura- 
tion before and after the December 27, 2004 burst. 



the azimuthal component of the magnetic field and in the 
electric field E, namely 




0.2 0.4 0.6 0.8 1 
Normalized magnetic fiux 



0.2 0.4 0.6 0.8 1 
Normalized magnetic fiux 



Fig. 6. The rescaled electric current distribution A/(i/'openi^Fo) 
and the rescaled distribution Aj4'/(i/'openi^^„), as functions of 
the rescaled magnetic flux i/'/i/'open in the open line region, for 
Q.Fo =1, 0.8 & 0.6 (from the lower curves up respectively). 
The upper curves (dotted) are the ones that correspond to the 
Michel split monopole expression. 



Q.Fo- Ea.l25limplies that stellar magnetic field estimates need to 
be revised upwards over the canonical value obtained when one 



r I Q.FrBp \ ^ Ql^Blrl (—] (—] (24)Compares eqs. 1 & 2. Note that when Qpo = 0, p,, = 0, 7 = 0, 
J \ I ync/ ync/ i.e. no currents flow through the magnetosphere, and therefore 



Here, the integration distance r extends to distances » ri^. Any 
evolution between the different solutions will require the re- 
lease (or buildup) of the corresponding energy difference (see 
discussion in the next section). 

We discovered that, as Q.fo varies from Q. to 0, the open 
field region decreases to a minimum value of about (Aopen ~ 1 .2 
(see fig.|5}. In the next section we will see that this numerical 
result might have interesting physical implications in under- 
standing the SGR phenomenon. 

Fig |6l shows the corresponding rescaled electric cur- 
rent distribution A/(i/'openf^Fo), and the rescaled distribution 
AA' /{tffopen^^o), (obtained numerically) as functions of the 
normalized magnetic flux i^/i^open- We see that indeed the elec- 
tric current distributions are very similar and proportional to 
Qfo. Let us now see how this result affects our estimation of 
stellar magnetic fields B,. As we mentioned in the introduc- 
tion, it is customary to estimate B, by equating the observed 
stellar spindown energy loss to the estimated electromagnetic 
spindown torque. As we show in the Appendix, 



-^em spindown 



Q 



Ji/<=0 



A((A)diA 



2 



open 



Bia'Q.For° /n, 
4c3 



(25) 



(in real units). In general, introduces one more free param- 
eter in the problem (see section 5). Let us here consider only 
the natural case rc ~ nc and discuss the physical significance of 



gnetosphe 

the star will not spin down. In most cases, Qfo ~ [80, 95]%Q 
(Romani, personal communication), and therefore, the correc- 
tion introduced in the stellar magnetic field estimate is in most 
cases practically insignificant. The correction is significant and 
should be taken into serious consideration for slow pulsars near 
the pulsar death-line, where y,(iAopen) ~ lO'^ Volts - V(ipopen) 
andVHfAopen) ^0(eq.[T3}. 

4. A 'coughing' magnetosphere 

The solutions presented in the previous section are all steady- 
state solutions characterized by one parameter, Qfo, which, as 
we argued, is determined by the particle acceleration gap mi- 
crophysics. Let's imagine first that charge carriers are freely 
available at the base of the magnetosphere. In that case, the gap 
is shorted out, and the magnetosphere is described by a steady- 
state solution with Qf „ ^ Q. (CKF). Let us imagine next that the 
supply of charge carriers is somehow suddenly depleted. The 
gap will suddenly grow, and the magnetosphere will quickly 
evolve towards a different steady-state solution with Q^o Q. 
We are now going to discuss how, in our opinion, the magne- 
tosphere may evolve from the one steady-state solution to the 
other. We will base our discussion on the particular example of 
SGR 1806-20, and its December 27,2004 burst. 

We will argue that, when the particle acceleration gap at the 
base of the magnetosphere suddenly grows, the magnetosphere 
will spontaneously evolve from a configuration with a larger 
open field line region and a larger poloidal electric current, to 
one with a smaller open field line region and a smaller poloidal 
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Fig. 7. Schematic magnetospheric evolution characterized as 
'magnetospheric coughing' (clockwise from upper left cor- 
ner). In the upper left corner is shown a steady-state solution 
with Qfo ~ Q.. Conditions at the base of the magnetosphere 
changed suddenly towards a different steady-state solution with 
Qfo 5^ O, and a spherical electromagnetic wave (shown with 
dotted line) sweeps through the open field line region at the 
speed of light. When the wave reaches the light cylinder region, 
reconnection allows the expulsion of the amount of magnetic 
flux required for the magnetosphere to evolve towards the new 
steady-state solution that corresponds to the new value of Qfo 
(upper right corner). The detached magnetic flux forms a 'plas- 
moid' fliat escapes in the equatorial region (lower right comer). 
The system reaches a final steady-state shown schematically in 
the lower left corner, and will remain there for as long as the 
physical conditions that sustained the new value of Q.fo at the 
base of the magnetosphere persist. 



electric current. One way to achieve this might be through 
north-south reconnection at the distance of the light cylinder. 
We expect a significant amount of magnetic flux (~ 5% i/'open) 
to 'snap' and move equatorially outward similarly to a solar 
coronal mass ejection (plasmoid). At the same time, the mag- 
netosphere will release the excess energy contained in the az- 
imuthal component of the magnetic field and in the elec- 
tric field E through a spherical electromagnetic wave sweeping 
through the open field region at the speed of Ught^. As is shown 
in eq. |23 the energy contained in that wave would grow with 
distance. We would like to characterize this dramatic evolution 
as 'magnetospheric coughing' (see fig. for a schematic de- 
scription). 

As long as the depletion of charges persists, the magneto- 
sphere will remain in the low Q.f„ state. The magnetosphere 
might return to a higher ilf,, state where angular momentum is 
removed more efficiently only if charge carriers become freely 
available again at the base of the magnetosphere. We specu- 
late that in such case, the magnetosphere will evolve through 



* In general, this will be a spherical Alfven wave moving outward 
at the Alfven speed. 



differential rotation between the star and the light cylinder re- 
gion, and therefore the evolution will be less dramatic than the 
magnetospheric coughing described above. 

In our example (see fig.|5|l, let us choose the solution with 
Q.Fo - 0.8 as the pre-burst solution. Our numerical analysis 
yielded 



Fo pre-burst 



0.8 



'Aopen pre-burst ~ 1 -23 



(26) 



(27) 



Based on our detailed axisymmetric ideal MHD model, and 
given the observed pre-burst spindown rate P = -8.5 x 
10"'^ Hz/sec, we obtain 



B, =5.2x 10''*G, 



(28) 



and a corresponding accelerating potential in the magneto- 
spheric gaps 



Vpre-burst = 5 X 10'° VoltS . 



(29) 



We know that, after the burst, the spindown rate was 2.7 times 
smaller JWoods et al. 20051 . This allows us to take 



^Fo post-burst — 0.31 , 
'Aopen post-burst — 1 -20 . 



V, 



post-burst 



1.6 X 10" Volts 



(30) 
(31) 
(32) 



We see that both before and after the burst, the accelerating po- 
tential is of the order of 10" Volts. Indeed, the magnetosphere 
is emitting pulsed radiation in both cases, only after the burst, 
pulsed radiation is observed to be weaker. We attribute this dif- 
ference to the smaller radiation cone (due to the smaller open 
field line region) which might thus avoid our line of sight. 

According to eq.|24] the energy difference between the two 
magnetospheres is of the order of 



^ Fo pre-burst 



" ^Fo post-burst)^* ''* ) ( H 



lO^M-) f-lergs. 



(33) 



According to|33l the energy contained in the spherical blast 
wave will be comparable to the apparent burst luminosities ob- 
served on earth (e.g. lYamazaki et al. 2005^ at distances r/ri^ ~ 
{r\clr»f. We would like to defer a more detailed discussion of 
the burst energetics to a future work. 

5. Conclusions 

In our present work we presented global numerical solutions 
of the generalized pulsar equation that describe the steady- 
state structure of axisymmetric rotating neutron star magneto- 
spheres. We have introduced two new parameters besides the 
neutron star angular velocity f2. 
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- Qfo, the angular velocity of rotation of open field lines. 
This quantity is related to the particle acceleration gaps at 
the base of the magnetosphere (the closer Q.fo is to Q., the 
smaller the gap), and is determined by gap microphysics 
outside the context of our present ideal MHD formulation. 

- rc, the maximum equatorial extent of the closed line region 
(see Appendix). We speculate that rc might be determined 
by inertial effects outside the context of our present ideal 
MHD formulation. 

Note that, in our global solutions, i/'open (the amount of open 
field lines) is determined self-consistently, and consequently it 
is not a free parameter (see, however, Goodwin et al. 2004 for a 
different point of view). Similarly to CKF, the poloidal electric 
current distribution that guarantees smoothness and continuity 
at the open field light cylinder is obtained itteratively, and an 
approximate analytic expression is given. Our results general- 
ize the solution presented in CKF: IGruzinov 20051 

We also obtained a generalized expression for the steady- 
state spindown magnetospheric energy losses (eq. \25l . which 
is different from the canonical one for a misallingned 
magnetic rotator Magnetospheres with different values of 
Qfo and/or rc contain different amounts of electric cur- 
rents, and therefore spin down differently. This changes 
slightly our estimates of stellar magnetic fields B, (see also 
[Harding, Contopoulos & Kazanas 1999] for a relevant discus- 
sion in the case of magnetar magnetic field estimates). More 
importantly, however, this might have serious implication in the 
calculation of the magnetic braking index n = Q.C1/Q?. One can 
easily check (ea. l25t that any functional dependence of Q.f and 
(Aopen different from the canonical one Q.f cc Q, and i/'open ^ 
will yield a braking index n 3 as obtained observationally 
(Contopoulos & Spitkovsky, in preparation). 

Finally, we argued that the magnetosphere may sponta- 
neously evolve between steady-state configurations character- 
ized by different values of Q.f„ and/or r^. The evolution from 
a high to low value of Q^o and/or low to high value of will 
result in the dramatic release of a significant amount of elec- 
tromagnetic field energy and magnetic flux. The return to the 
former configuration will be less dramatic, since it will require 
the buildup of the corresponding electromagnetic field energy 
difference. Our results might be relevant in understanding the 
SGR burst phenomenon. 
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Appendix A: Pulsar spindown estimates 

When a neutron star is not surrounded by vacuum, the rotating 
charged relativistic Goldreich- Julian-type magnetosphere is threaded 
by poloidal and toroidal electric currents. We will consider only the 
axisymmetric case for simplicity. Two large scale poloidal electric cur- 
rent circuits (north & south) are generated. These flow only along open 
field lines, and close along the surface of the star at the two polar caps 



where they generate electromagnetic torques antiparallel to the angu- 
lar momentum of the neutron star 



1 



-rBJdSdr 



(A.l) 



through any stellar cross section d5 threaded by poloidal electric cur- 
rent density J. One can easily check that the stellar kinetic energy loss 
through the above torques is given by 



■^em spindown 



: Q. 



^=1 



A(l/r)dl/r - -Q.Q.p,^l 



(A.2) 



(our expression accounts for the two hemispheres, north & south). We 
made use of the numerical result i^opcn ~ 1-23. At the same time, the 
magnetosphere radiates electromagnetic energy 



An 



rEpB^AS 



(A.3) 



through any cross section AS in the region of open field lines. One 
can easily check that the total electromagnetic energy loss through the 
above Poynting flux is given by 



\p=0 



(A.4) 



flf is in general smaller than f2, and therefore, L^^ is in general less 
than spindown- The difference between the two is consumed in the 
particle acceleration gaps that develop along open field lines, namely 



^paiticles ^em spindown ^em ' 



A{ilj){Q. - np(if/))dif/ 



(A.5) 



The above expressions are normalized to the Goldreich- Julian value 



4c3 



(A.6) 



Appendix B: Alternative magnetospheric solutions 

In solving eg. 1101 we have all along argued that nature will choose 
the most natural solution, namely the one with the maximum extent 
of the 'dead zone'. A competing to the above scenario might be one 
where the extent of the 'dead zone' is a free parameter (Goodwin et 
al. 2004). Since Q.f introduces one more free parameter in the prob- 
lem, we will consider only one representative case with Q.fo - 0.8. 

can be solved numerically as described before. In this sce- 
nario, solutions with a smaller 'dead zone' are also more efficient in 
removing angular momentum from the spinning star fea. lA.2> . As an 
example, we take ~ 0.59 and obtain i/'open = 2.03 (fig. lB.H . This so- 
lution may evolve rapidly through reconnection towards the solution 
shown in fig. |2| with ~ I, i/'open = 1.23, and thus yield a spindown 
rate 2.7 times lower, releasing at the same time a significant amount 
of magnetic field energy. Note that the system is an efficient radiator 
through particle acceleration processes both before and after the burst 

(ilfo < 1). 
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